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AEROELASTICITY OF WING AND WING-BODY 
CONFIGURATIONS ON PARALLEL COMPUTERS 


I. INTRODUCTION 

There is a continuous effort to improve the performance of subsonic transport aircraft. 1 
One attempt is to improve the fuel efficiency by extending the flight regime to high sub- 
transonic Mach numbers to increase lift-to-drag ratios and flight speeds. To avoid the 
high drag associated with strong shock waves, these advanced transports require modern 
wing sections such as supercritical wings that delay the shock wave formation. Early 
experiments have shown that these advanced wings experience an undesirable reduction 
in the flutter speed at the transonic regime. Such a phenomenon commonly known as 
‘transonic dip’, is more pronounced for wings with supercritical airfoils. 2 Furthermore, for 
accurate prediction of the flutter characteristics, it is necessary to model viscous flows 
using the Navier- Stokes equations. 

To date, advanced wing calculations have been limited to steady and unsteady compu- 
tations on rigid wings. However, it is necessary to account for the structural flexibility to 
accurately compute its aeroelastic characteristics. The aeroelastic deformation resulting 
from this flexibility can significantly change the nature of the flow. Strong interactions 
between the flow and structures can lead to sustained aeroelastic oscillations for swept 
wings. 3 Also, it is necessary to include the flexibility for proper correlations of computed 
data with experiments, particularly with those obtained from flight tests. To compute 
the flows accurately, it is necessary to include both aerodynamic and structural effects of 
the body. Recent efforts have been made to include the flexibility effects for wing-body 
configurations. 4 

Furthermore, in recent years, significant advances have been made for parallel com- 
puters in both hardware and software. Now parallel computers have become viable tools in 
computational mechanics. In order to exploit the new architecture of parallel computers, 
the computer code, ENSAERO, 5 which computes the unsteady aerodynamics and aeroe- 
lasticity of aircraft by using the thin layer Navier-Stokes equations, has been parallelized on 
the Intel iPSC/860 parallel computer. The parallel version of the code has demonstrated 
its capability to compute aeroelastic responses by concurrently integrating the Euler equa- 
tions and the structural equations of motion with aeroelastically deforming grids. 6 The 
code can also model unsteady viscous flows using time-accurate, finite-difference schemes 
based on the Beam-Warming algorithm with the Baldwin-Lomax turbulence model. 

The objective of this research is to develop computationally efficient methods for solv- 
ing aeroelasticity problems on parallel computers. Both uncoupled and coupled methods 
are studied in this research. For the uncoupled approach, the conventional U-g method is 
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used to determine the flutter boundary. The generalized aerodynamic forces required are 
obtained by the pulse transfer-function analysis method. For the coupled approach, the 
fluid-structure interaction is obtained by directly coupling finite difference Euler/Navier- 
Stokes equations for fluids and finite element dynamics equations for structures. This 
capability will significantly impact many aerospace projects of national importance such 
as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin be- 
comes very critical at the transonic region. This research effort will have direct impact on 
the High Performance Computing and Communication (HPCC) Program of NASA in the 
area of parallel computing. 

II. PREVIOUS STATUS 

A multidisciplinary code for computing unsteady flows and aeroelastic responses of 
aerospace vehicles, ENSAERO, has been developed on serial supercomputers at the Ap- 
plied Computational Aerodynamics Branch of the NASA Ames Research Center . 7 This 
multidisciplinary code computes unsteady aerodynamic responses of aircraft using the 
Euler /Navier- Stokes equations. An aeroelastic shape-conforming moving grid is used to 
include the effect of structural deformations on unsteady flows. This code is designed 
in a modular fashion to adopt several different numerical schemes suitable for accurate 
aeroelastic computations. The basic coding of ENSAERO can accommodate zonal grid 
techniques for efficient modeling of full aircraft. 

An early version of ENSAERO 8 has been successfully applied in computing aeroelastic 
responses of a rectangular wing by using the Euler equations for fluids and the modal 
equations for structures. The result demonstrates that the code can accurately predict 
the flutter dynamic pressure of a rectangular wing. The code was extended to compute 
aeroelastic responses using the Navier-Stokes equations for fluids . 9 Later, it was updated by 
utilizing an upwind algorithm, and the code has been applied to fighter wings undergoing 
unsteady motions 10 ’ 11 at moderately large angles of attack. This code also has a capability 
of modeling moving control surfaces . 12 Furthermore, ENSAERO has demonstrated the 
capability to simulate transonic flows on wing-body configurations using the Navier-Stokes 
equations . 13 

In the past, the modal equations were used to model structures for the purpose of 
aeroelastic analysis. For simple geometries such as clean wings, the modal approach can 
predict accurate response results. However, the modal approach may be less accurate for 
complex structures such as wing-body configurations. In order to accurately represent 
aeroelastic responses of general wing-body configurations, the modal equations should be 
replaced with the finite element equations. Recently, a typical wing-body configuration 
has been used to demonstrate aeroelastic responses at transonic Mach numbers using 
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the Navier-Stokes equations for fluids and the finite element equations for structures. 14 
Simple one-dimensional beam elements are used to model the wing-body structures. Each 
node has three degrees of freedom (DOF) corresponding to transverse displacement and 
to transverse and torsional rotations, respectively. 

Recently, a version of ENSAERO 15 that uses the Euler equations for fluids and the 
modal equations for structures has been parallelized on the Intel iPSC/860 at Ames. The 
Intel iPSC/860 is a distributed-memory, multiple-instruction, multiple-data (MIMD) com- 
puter with 128 processors. In this parallel implementation, a domain decomposition ap- 
proach is used in which the fluid equations and the structural equations are modeled in 
separate computational domains. Each domain is mapped individually onto a group of 
processors, referred to as a cube on the Intel iPSC/860. However, because of the coupling 
between the disciplines, there is a need to exchange data, such as pressures and structural 
deformations at interfaces. This exchange between the fluid and structural domains is 
accomplished through an intercube communication mechanism, 16 which enables different 
processors in each cube to communicate directly. 

III. CURRENT WORK 

1. FINITE ELEMENT MODELING OF STRUCTURES 

The finite element representation of structures generally provides more accurate mod- 
eling of structures than the modal representation does in aeroelastic computations. For 
the aeroelastic computations of wing and wing-body structures, plate and shell models are 
used as reported in References 4 and 17. The parallel implementation of the structural 
domain is mainly discussed in Ref. 17. Although the current implementation is only us- 
ing ANS4 plate/ shell elements to model structures, it is possible to use different types of 
elements. 

2. FLUID-STRUCTURAL INTERFACE 

In aeroelastic analysis, it is necessary to represent the equivalent aerodynamic loads at 
the structural nodal points and to represent the deformed structural configurations at the 
aerodynamic grid points. In the present domain decomposition approach, coupling between 
the fluid and structural domains is achieved by interfacing the boundary data, such as 
aerodynamic pressures and structural deflections, at each time step. An analytical moving- 
grid technique has been successfully used to deform the aerodynamic grid according to the 
structural deflections at the end of every time-step. 7,8,14 There are different approaches 
for obtaining the external load vector, depending on the equations used for the structural 
dynamic analysis. 
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In order to replace modal equations with finite element structural dynamics equations 
for fluid-structural interaction problems, a new fluid-structural interface, similar to the 
modal matrix used for modal equations, should be developed. Several numerical proce- 
dures have been developed for exchanging the necessary information between the fluid and 
structural domains. 18-20 In this research, two different types of fluid-structural interfaces 
were studied and compared as shown in Ref. 4. 

The current implementation of the fluid-structural interface on the Intel iPSC/860 
is based on direct node-to-node communication. Thus each of precessors assigned to the 
fluid domain can communicate with any processors of the structural domain. A typical 
communication patera is generally irregular during the aeroelastic computation of a High 
Speed Civil Transport (HSCT) model on the Intel iPSC/860. More details are found in 
Ref. 6. 

3. PARALLEL INTEGRATION 

In a serial computer, the integration of both fluid and structural equations is performed 
one after the other in a sequential nature. When implementing the integration scheme on 
parallel computers, all processors can be used to solve the fluid and structural equations 
sequentially. But this approach requires more memory per processor and two disciplines 
have to be implemented in a single program. As a result, modularity of each algorithm 
for individual disciplines will have to be sacrificed to a significant degree. In addition, this 
approach will be less efficient as increasing the number of processors because the problem 
is not linearly scaled. 

However, while keeping modularity of each discipline, computations can be done more 
efficiently on MIMD parallel computers by executing the integration of both fluid and 
structural equations concurrently. In the proposed parallel integration scheme, both do- 
mains start computations independently and one of the solvers waits until the other finishes 
its calculation. Then they exchange the required data with each other for the next time 
step. By doing so, the parallel integration can reduce the idle time since only one cube 
(the fastest) will have to wait. This integration scheme exploits the parallelism offered 
by the domain decomposition approach to solve the coupled fluid-structural interaction 
problems. More details are found in Ref. 6 and 17. 

4, STATIC AEROELASTICITY 

It is of interest to find static deflections of aerospace vehicles subjected aerodynamic 
loads. In the past, static equilibrium equations of structures are used in ENSAERO to find 
the static deflections in conjunction with the aerodynamic forces computed by spacially 
variable time method starting from the steady state flow solutions as initial conditions. 
This method, however, shows highly oscillatory behavior in the solutions. Sometimes this 
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approach becomes highly unstable and difficult to obtain solutions. In order to overcome 
these problems, dynamic equilibrium equations with high dampimg are suggested to re- 
place the static equilibrium equations of structures. 21 The matrix from of the dynamic 
equilibrium equations of motion is 


[M]{q] + [G]{q} + [K]{q} = { Z} 


(1) 


where [M], [(?], and [A'] are the mass, damping and stiffness matrices, respectively. {Zj is 
the aerodynamic force vector. By adding the inertia and damping terms, the variation of 
displacements becomes smooth and the final solution approaches to the static equilibrium 
positon of structures. This approach can start from the free stream conditions and obtain 
the steady state flow solution at the deformed configuration in addition to the static 
deflections. Furthermore, as pointed out in Ref. 21, this approach is quite stable and 
efficient compared to the previous approach. This method has been implemented in the 
current parallel version of ENSAERO. 

5. FORMULATION FOR FLUTTER BOUNDARY PREDICTION 


The direct coupling approach requires unprecedented amount of computation time to 
determine the flutter boundary in the past. Thus uncoupled approaches have been widely 
used to determine the flutter boundary. One of the uncoupled approaches is the U-g 
method, which computes the flutter speed by solving an eigenvalue problem. Following 
the Ref. 22, the eigenvalue equations can be derived as 


[K]^[{M} + C[Q]]{q} = \{q} 

,3 


c = 


A = 


pc 

2k 2 

(1 + ig) 

u > 2 


complex eigenvalue 


( 2 ) 


where the reduced frequency, k = uc/U ) c is the reference chord, and U is the flight speed. 
The generalized aerodynamic forces are obtained as 


Qij = \J J h(x,y)iAC p (x,y)jdxdy 

h = modal displacement 


( 3 ) 


where i and j represent mode shapes. The coefficient Qij represents the force acting in 
the i th mode due to pressure generated by the unsteady motion of the j th mode and is 
dependent on the reduced frequency. For each reduced frequency chosen, the generalized 
aerodynamic force matrix [Q] is obtained by solving the unsteady flow equations and then, 
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in turn, the complex eigenvalue problem can be solved to obtain the damping coefficient, 
g. The damping coefficient, g, is used to determine the flutter point. 

One way to compute GAFs is the time integration (TI) approach. The TI method 
performs several cycles of a forced harmonic oscillation and uses the last cycle of oscillation 
to determine the first harmonic component of GAFs. Multiple unsteady computations are 
required at various reduced frequencies to generate GAFs for each structural mode used 
in the flutter analysis. Typically unsteady data at 3 to 4 given frequencies are required to 
make a reasonable prediction of the flutter point at a given flight speed when using the U-g 
method. Interpolation techniques are usually used to obtain the generalized aerodynamic 
forces at various reduced frequencies to accurately determine the flutter point. Although 
there have been great advances in conventional and parallel supercomputers, it is still 
formidable task to determine the flutter boundary using the time integration approach. 

Repetitive computations of GAFs for each frequency can be avoided by using the 
indicial approach 23 or pulse transfer-function analysis (PTFA) method. 24,25,26 In these 
approaches, GAFs for all frequencies for a given structural mode can be extracted from 
a single unsteady response computation. These approaches require an assumption that 
the unsteady flow can be linearized about a nonlinear steady flow. Such an assumption is 
valid for small perturbation. It is noted that classical flutter starts as a small perturbation 
phenomenon and, therefore, it is considered to be appropriate to use these approaches to 
make preliminary prediction of the flutter boundary. 

The indicial method computes GAFs from the Fourier transform of the indicial re- 
sponse. That is, a step change in the structural mode is used as the initial condition and 
the unsteady computation is continued until the transient had decayed. The PTFA method 
is a variation of the indicial method. The Fourier transform of the resulting force response 
is divided by the transform of the generalized displacement of the given structural mode 
to obtain GAFs. It uses a forced modal motion represented by a smoothly varying pulse 
function instead of a step function. Both approaches can significantly reduce the compu- 
tation time to generate GAFs compared to the time integration method. In this work, 
the generalized aerodynamic force matrix is computed by using both time integration and 
pulse transfer-function methods for comparison. However, the PTFA method is used to 
determine the flutter boundary. 

6. PARALLELIZATION OF FLUTTER BOUNDARY PREDICTION 

For the U-g method used to determine the flutter boundary, the generalized aerody- 
namic force matrix [Q] is required as functions of the reduced frequency. Although the 
U-g method is an efficient method to determine the flutter boundary, it would become a 
CPU-intensive procedure when generalized aerodynamic forces (GAFs) are computed by 
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solving the Euler or Navier-Stokes unsteady flow equations instead of simplified unsteady 
aerodynamic equations. So the parallelization is concentrated on computing GAFs. 

The parallelization of the flutter boundary prediction procedure can be achieved in 
three different ways. The first is the concurrent execution of unsteady computations as- 
sociated with different structural mode shapes. This is represented as a box with broken 
line in Fig. 1. This computation can be accomplished for a fluid grid which fits on a 
single processor node. Upon completing unsteady computations, the generalized aerody- 
namic force matrix can be assembled with the minimal communication among processors 
assigned to a given flight condition. Then the flutter point can be determined from the 
eigenvalue analysis. This can be repeated for several flight Mach numbers to obtain the 
flutter boundary. This is a coarse grain parallelization since parallelization is accomplished 
only for the given number of modal motions. 

The second is the concurrent multiple execution of the first case with different initial 
steady-state conditions, e.g. different Mach numbers. Each column in Fig. 1 represents the 
first approach and the multiple columns represent the execution of the multiple unsteady 
flow analyses in parallel. In this way, the flutter boundary for the given flight regime 
can be determined in a single run. Upon completion of unsteady modal motion, the pulse 
transfer-function analysis yields GAFs as functions of the reduced frequency. Flutter points 
can be determined independently by solving the eigenvalue equations for various reduced 
frequencies after assembling the GAFs. In order to assemble GAFs, communication for 
unsteady data should be limited only among the processors with the same initial condition. 
This is accomplished by using the MPIRUN, 27 a utility developed at NASA Ames Research 
Center. The MPIRUN enables to flexibly define a group of processors and communicate 
data within a group or between groups. The MPIRUN is based on the message passing 
interface (MPI) standard. 28 

However, since each unsteady computation is still quite CPU-intensive task, it is not 
easy to complete the whole simulation in a single run. Therefore, in order to cut down 
the total run time, it is necessary to parallelize each unsteady flow computation. This 
parallelization is also necessary for bigger problems where the fluid grid can not fit on a 
single processor node. This involves the parallelization of the flow solver and the moving 
grid scheme. This is referred to a fine grain parallelization. It requires another level of 
communication so that each modal motion can be subdivided and distributed onto M 
number of processors as shown in Fig. 1. Thus, during the solution stage of unsteady 
flow equations, communication should be limited among the processors with the same 
forced modal motion and the initial condition. More details for parallelization of the flow 
solver on the Intel iPSC/860 computer can be found in Reference 29. The allocation of 
processor nodes is determined by the number of vibration modes selected, the size of the 
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computational grid, and the number of cases for different flight speeds. By combining the 
above three approaches together, the procedure can determine the flutter boundary for a 
given wing with significant reduction in computational time. 

7. COMPUTATIONAL RESULTS 


a. Dynamic Aeroelasticity 

Dynamic aeroelastic computations for the NASA Langley Clipped Delta Wing were 
performed first on a Cray Y-MP serial computer and then on the Intel iPSC/860 MIMD 
parallel computer. Comparison between the modal and finite element analyses for struc- 
tures and accuracy of the interfaces considered in this research are given in Ref. 4. 

Computational performance results are found in Ref. 17. Including the data exchange 
between fluid and structural domains, the current aeroelastically deforming grid scheme 
requires about 12 percent of the computational time per each integration step. 

Application of the procedure to wing-body configurations can be found in Ref. 4 and 
6. Dynamic aeroelastic computations are done on flexible wing and body structures such 
as a simple wing-body and a Boeing 1807 HSCT model, the results are given in Ref. 6. 

b. Static Aeroelasticity 

For demonstration purpose of static aeroelasticity, the LANN wing was selected since 
this wing represent the ASCT type wing configuration. The LANN model has an aspect 
ratio of 7.92, a taper ratio of 0.40 and a leading-edge sweep angle of 27.5 degrees. The 
airfoil sections axe supercritical with constant thickness-to-chord ratio of 12 percent. This 
model was tested in the transonic wind tunnel (HST), the Netherlands. Details of steady 
flow measurements are given in Ref. 30. 

Computations are made using the built-in C-H grid topology available in the EN- 
SAERO code. For the flow computations, a grid of size 151 x 35 x 35 ( 151 points in 
streamwise direction, 35 points in both spanwise and surface-normal directions) is selected 
as it was found to be appropriate for this problem in the early work. 31 Using this grid, 
a steady computation is made at M = 0.87, Re = 7.6 x 10 6 degrees and AoA = 3.0 . 
As shown in Fig. 2, it is also noted that, the Euler equations did not properly represent 
the flow over the LANN wing. This figure shows that the Navier-Sokes equations with 
the Baldwin-Lomax turbulence model obtained the result in good agreement with the 
experiment. It is of interest to investigate the effect of viscosity on the static deflection. 

Figure 3 shows surface pressure distributions obtained using the Navier-Stokes equa- 
tions at the same free stream condition. The effect of viscosity is shown at Fig. 4. For 
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this particular problem, the leading edge deflection of the wing computed by using the 
Navier-Stokes equations are smaller than that obtained by using the Euler equations. 

Figures 5 and 6 shows the effect of flexibility on pressure distributions. Both figures 
show quite different trends when the wing is allowed to be flexible. For the Euler com- 
putations, the shock on the upper surface moves forward to the leading from the rigid 
solution whereas, for the Navier-Stokes computations, the shock position moved backward 
from the rigid solution. Figure 7 shows the surface pressure distributions at the deformed 
configurations . 

c. Parallel Flutter Boundary Prediction 

In order to validate the present procedure and compare the pulse transfer-function 
analysis (PTFA) method with the time integration (TI) method, computations were made 
for a unswept rectangular wing of aspect ratio 5 with a 6% thick circular-arc airfoil. 31 The 
model is 11.50 in. in the spanwise direction and in the 4.56 in. chordwise direction. The 
transonic flutter characteristics of this wing are available from wind tunnel tests for various 
flow parameters. For this computation, the flow field is solved using the Euler equations 
with a C-H grid of size 151 x 25 x 30. 

Generalized Aerodynamic Forces 

To compare the PTFA method with the TI method, the unswept rectangular wing was 
selected. The first three structural modes are used in the flutter analysis. The maximum 
amplitude of the generalized displacements is controlled to be small enough so that the 
response of the aerodynamic forces is linear. All computations were made on the Cray 
C90 computer for the purpose of comparison between the TI and PTFA methods. For 
the TI method, the unsteady data is obtained by integrating the unsteady flow equations 
when the wing is oscillating in harmonic motion for a given structural mode. The TI 
method required three cycles with 3600 time steps per cycle for a complete calculation. 
This computation required about 2 CPU hours on a single Cray C-90 processor. Using the 
PTFA method, the unsteady data are obtained by perturbing the wing with a smoothly 
varying pulse input for a given structural mode. The total CPU time required for a single 
PTFA run is about 1.5 hours on a single Cray C-90 processor. It should be noted that the 
unsteady computation has to be continued until the transient has decayed. 

Figures 8 and 9 show the lift coefficient history obtained by using the TI and PTFA 
methods, respectively. Both results are obtained from unsteady computations by using the 
Euler equations at free stream Mach 0.715 and zero angle of attack. For the TI method, 
since a harmonic motion is enforced, the residual shows that the transient dies out and 
the lift coefficient show periodicity in response as shown in Fig. 8. Since a pulse motion 
is enforced for the PTFA method, there is no periodicity in response as shown in Fig. 
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9. However, the lift coefficient response from the PTFA method shows that the transient 
due to the pulse motion has decayed out and the solution converged to the steady-state 
condition. 

The comparison of the GAFs obtained using the TI and PTFA methods is presented 
in Fig. 10. The results are obtained for the unswept rectangular wing at the free stream 
Mach 0.715 and zero angle of attack by using the Euler flow equations. Both results show 
in good agreement. It was observed in the past that the results from the PTFA method 
became deviated from the results obtained by the TI method as increasing the reduced 
frequency. 33 However, as shown in Fig. 10, it is considered that the PTFA is applicable to 
generate GAFs within the range of the reduced frequency required for the flutter prediction. 

Flutter Determination 

The present flutter prediction procedure has been validated by comparing the com- 
putational results with the experiment. The comparison between the results from the 
computation and experiment is made in the Fig. 11. The computed flutter frequencies 
axe compared with those obtained in the experiment. The computational results show 
similar trend as compared to the experiment at the lower Mach number range although 
the computed result predicts the transonic dip earlier than the experiment. The difference 
in the prediction of the transonic dip can be attributed to the possible viscous effect since 
the Euler equations are used in the computation. 

IV. SUMMARY AND CONCLUSIONS 

In this work, procedures to compute aeroelastic responses of wing and wing-body 
structures have been parallelized on MIMD parallel computers using uncoupled and cou- 
pled approaches. For the coupled approach, dynamic aeroelastic responses are obtained 
by direct coupling of the finite difference Euler/Navier Stokes equations for fluids and the 
finite element dynamics equations for structures for wing and wing-body structures. The 
procedure is based on a domain decomposition approach which enables algorithms for the 
fluid and structural disciplines to be developed and maintained independently. 

A parallel integration procedure is developed to solve fluid and structural equations 
together. The parallel integration scheme enables the combination of advanced CFD and 
CSD technologies with minimal increase in computational time per integration step while 
keeping modularity of each discipline. The time per integration step is solely determined 
by the domain that requires most computational time on the iPSC/860. This parallel 
integration is one of the advantages of using MIMD computers for multidisciplinary anal- 
ysis. The procedure developed in this research will provide an efficient tool for solving 
aeroelastic problems of complete aerospace vehicle configurations on MIMD computers. 
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For the uncoupled approach, the conventional U-g method is used to determine the 
flutter boundary. Current implementation of the proposed procedure for the flutter predic- 
tion using the U-g method is limited to coarse grain parallelization. Thus the flow solver 
is running on a single node and several flow conditions and oscillating modes are executed 
simultaneously. 

Based on this work the following conclusions can be made. 

1. It is feasible to directly couple the finite difference flow equations and finite 
element structural equations to obtain accurate results, though each discipline 
is solved in a separate computational domain. This domain decomposition ap- 
proach takes advantage of efficient methods developed for each individual dis- 
cipline. This approach can be extended to include more other disciplines, e.g., 
controls, optimizations, etc. 

2. The use of pulse transfer-function method for generalized aerodynamic forces 
(GAFs) is suitable for the problems solved in this work. This method is very 
efficient to obtain GAFs as a function of the reduced frequency. 

3. The parallel procedure for the prediction of the flutter boundary can obtain the 
flutter boundary in efficient manner on MIMD parallel computers. 
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Steady Cp 



Figure 1 . Comparison of sectional pressure coefficients among the 

computations and experiment at the freestream condition of M 
0.871 , angle of attack 3.0 degrees, Reynolds number 7.472 
million. (Code : ENSAERO, Grid : 151x35x35) 
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Figure 2. Surface pressure distributions obtained using Navier- Stokes 
equations at the freestream condition of M = 0.871 , angle of 
attack 3.0 degrees, Reynolds number 7.472 million. 

(Code : ENSAERO, Grid : 151x35x35) 
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Figure 3. Effect of viscosity on static aeroelastic responses computed by 
directly coupling the Navier-Stokes flow equations wltn finite 
element structural dynamics equations on the Intel iPSC/860 
computer (36 processors). 

(CFD Grid : 151x35x35, FE model : 128 plate elements) 
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Figure 4. Effect of flexibility on pressure distributions of steady state 
solutions obtained by using the Euler equations. 

(CFD Grid : 151x35x35, FE model : 128 plate elements) 



Figure 5. Effect of flexibility on pressure distributions of steady state 
solutions obtained by using the Navier-Stokes equations. 
(CFD Grid : 151x35x35, FE model : 128 plate elements) 




Euler computation 



Navier-Stokes computation 


Figure 6. Surface pressure distributions on deformed configurations 
obtained at the freestream condition of M = 0.871 , angle of 
attack 3.0 degrees, Reynolds number 7.472 million. 

(Code : ENSAERO, Grid : 151x35x35) 
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Figure 7. A schematic diagram for the procedure of flutter boundary prediction using the U-g method 




















Nondimensional time 

Figure 8. Time history of the total lift coefficient of the unswept rectangular 
wing with 6% thick circular-arc airfoil under a harmonic motion. 
(M = 0.71 5, angle of attatck = 0.0, Euler) 



Nondimensional time 

Figure 9. Time history of the total lift coefficient of the unswept rectangular 
wing with 6% thick circular-arc airfoil under a smoothly varying 
pulse motion. (M = 0.715, angle of attatck = 0.0, Euler) 
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Reduced Frequency 

Figure 10. Comparison of generalized aerodynamic forces obtained by using the time integration 
and pulse transfer-function analysis methods. (Code : ENSAERO, Grid : 151x35x35) 
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Figure 1 1 . Comparison of the computed flutter frequency of an rectangular 
unswept wing with the experiment. 











